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ABSTRACT 

The recent Hubble Space Telescope near-infrared imaging with the Wide-Field Camera 
#3 (WFC 3) of the GOODS-South field in the CANDELS program covering nearly 
lOOarcmin 2 , along with already existing Advanced Camera for Surveys optical data, 
makes possible the search for bright galaxy candidates at redshift z pa 7 — 9 using 
the Lyman-break technique. We present the first analysis of z'-drop z w 7 candidate 
galaxies in this area, finding 19 objects. We also analyse y-drops at z w 8, trebling the 
number of bright (Hab < 27 mag) V-drops from our previous work, and compare our 
results with those of other groups based on the same data. The bright high redshift 
galaxy candidates we find serve to better constrain the bright end of the luminosity 
function at those redshift, and may also be more amenable to spectroscopic confirma- 
tion than the fainter ones presented in various previous work on the smaller fields (the 
Hubble Ultra Deep Field and the WFC 3 Early Release Science observations). We also 
look at the agreement with previous luminosity functions derived from WFC 3 drop- 
out counts, finding a generally good agreement, except for the luminosity function of 
Yan et al. (2010) at z pa 8, which is strongly ruled out. 

Key words: galaxies: evolution galaxies: formation galaxies: starburst galaxies: 
high- redshift ultraviolet: galaxies 



1 INTRODUCTION 

Thanks to the installation of Wide Field Camera 3 (WFC3) 
on the Hubble Space Telescope (HST) in Summer 2009, the 
search for star-forming galaxies at redshifts z > 7 with the 
Lyman break technique (see Section 3) has become possi- 
ble with the infrared channel of the Wide-Field Camera #3 
(WFC 3) and led to the discovery of several galaxy candi- 
dates at z ~ 7 — 10. From these candidates we can deter- 
mine the rest frame UV luminosity function (LF) at these 
redshifts (Bunker et al. 2010, Wilkins et al. 2011a, Loren- 
zoni et al. 2011, Bouwens et al. 2011), an important tool 
in understanding the star formation history of the Universe, 
and also crucial to addressing the role of star-forming galax- 
ies in reionization. These works show a broad agreement on 
the clear LF evolution from z — 6 (and below) to z — 7 
with the characteristic luminosity L* fainter at higher red- 
shifts, and suggest further evolution at even higher redshifts 

* E-mail: silvio.lorenzoni@astro.ox.ac.uk 



(z ~ 8 — 10), although based on fewer candidates. The wealth 
of WFC 3 data on the Great Observatories Origins Deep 
Survey South (GOODS-S) area recently obtained by The 
Cosmic Assembly Near-infrared Deep Extragalactic Legacy 
Survey (CANDELS, Grogin et al. 2011, Koekemoer et al. 
2011), covering an area twice as large as the area surveyed 
in our previous papers, allows us to put better constraints on 
the bright end of the UVLF at z ~ 7— 9. The larger field now 
available also allows the identification of brighter sources, 
which may be more amenable to spectroscopic follow-up. 
The fact that these new WFC 3 images coincide with exist- 
ing deep Adavanced Camera for Surveys (ACS) optical im- 
ages is critical in rejecting potential interlopers - the ACS 
filters lie below the Lyman limit and hence any detection 
at short wavelength will reject low redshift contaminants. 
This is a luxury not afforded to recent pure-parallel surveys 
for high-redshift drop-outs like the Hubble Infrared Pure 
Parallel Imaging Extragalactic Survey (HIPPIES, Yan et al. 
2011) and the Brightest of Reionizing Galaxies (BoRG) sur- 
vey (Trenti et al. 2011 & Bradley et al. 2012). In this pa- 
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per we present for the first time a list of z'-drops at z ~ 7 
drawn from the large CANDELS field of GOODS-S. We also 
present our selection of z ~ 8 Y-drops in this field, and com- 
pare this with recent independent analyses of CANDELS 
Y-drops in GOODS-S by Oesch et al. (2012) and Yan et al. 
(2012). 

This paper is organised as follows: in Section 2 we out- 
line the HST observations with WFC 3 and the data reduc- 
tion, and in Section 3 we describe our colour selection to re- 
cover high-redshift Lyman break galaxies, and compare our 
sample with those from other studies. In Section 4 with dis- 
cuss the UV luminosity function derived from the new data. 
Our conclusions are presented in Section 5. Throughout, we 
adopt the standard concordance cosmology of Qm = 0.3, 
Qa = 0.7 and use Ho = 70kms _1 Mpc -1 . All magnitudes 
are on the AB system (Oke & Gunn 1983). 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Observations 

In this paper we analyse images from WFC 3 on HST taken 
in the F105W, F125W and F160W filters, corresponding ap- 
proximately to the near-infrared Y-, J- and H -bands. The 
data come from the HST programs GO-12060, GO-12061 
and GO-12062 in the CANDELS program (P.I. S. Faber, 
see Grogin et al. 2011; Koekemoer et al. 2011), covering the 
areas of the GOODS-S field (Giavalisco et al. 2004) not cov- 
ered by the Early Release Science (ERS) program GO /DD- 
11359 (P.I. R. O'Connell, see Wilkins et al. 2010). The area 
is divided into a 'deep' field, measuring ~ 63 arcmin 2 with 

3 orbits in each Yiosw, Ji25w and //i60w filters, and a 'wide' 
field with one orbit per filter over an area of ~ 33 arcmin 2 
(the areas quoted refer to the deepest area where the cover- 
age has the maximum number of overlapping frames). Ex- 
tensive ACS imaging has been carried in these areas in previ- 
ous years (Giavalisco et al. 2004, Beckwith et al. 2006) in the 
6 (F425W), v (F606W), i (F814W) and z' (F850LP) filters, 
allowing us to confidently use the Lyman-break technique 
to select likely high redshift star forming galaxies. 

The infrared channel of WFC 3 was used, which is a 
Teledyne 1014 x 1014 pixel HgCdTe detector (a 10-pixel 
strip on the edge is not illuminated by sky and used for 
pedestal estimation), with a field of view of 123" x 136". 
The data were taken in "MULTIACCUM" mode using 
SPARSAMPLE100, which non-destructively reads the array 
every 100 seconds. These repeated non-destructive reads of 
the infrared array allow gradient-fitting to obtain the count 
rate ("sampling up the ramp") and the flagging and rejec- 
tion of cosmic ray strikes. In Table [1] we list the exposure 
time for both the 'deep' and 'wide' fields for each spectral 
band. 



2.2 Data Reduction 

Data reduction is performed as described in our previous 
papers (Lorenzoni et al. 2011, Wilkins et al. 2011a). We 
used the IRAF.STSDAS pipeline calwf c3 to calculate the 
count rate and reject cosmic rays, then MULTIDRIZZLE 
(Koekemoer et al. 2002) to combine exposures taking ac- 
count of the geometric distortions and mapping on to an 



output pixel size of 0'.'06 from an original 0'.'l3pix _1 , which 
corresponds to a 2 x 2 block-averaging of the GOODSv2.0 
ACS drizzled images in b-, v-, i- and z'-bands. We used a 
MULTIDRIZZLE pixel fraction of 0.8 for the 'deep' area 
and 1.0 for the 'wide' area to recover some of the under- 
sampling. We used our own reduction of all the WFC3 data 
for the CANDELS GOODS-S 'wide' area and of the Y-band 
data of the 'deep' region. For the J- and //-bands covering 
the 'deep', we used the reduced single epoch images made 
available by the CANDELS teanQ and co-added these to- 
gether with inverse-variance weighting (i.e. weighting each 
pixel by its exposure time). 

For WFC3, we use the zeropoints reported on 
|http : //www, stsci . edu/hst/wf c3/phot _zp_lbn| last up- 
dated in January 2011, where the zeropoints are 26.27, 26.25 
& 25.96 for F105W, F125W & F160W. 

We perform photometry using fixed apertures of 0'.'6 
diameter, and introduce an aperture correction to account 
for the flux falling outside of the aperture. This correction 
was determined to be ~ 0.2 — 0.25 mag in WFC3 from pho- 
tometry with larger apertures on bright but unsaturated 
point sources. For the ACS images, the better resolution 
and finer pixel sampling require a smaller aperture correc- 
tion of « 0.1 mag. All the magnitudes reported in this pa- 
per have been corrected to approximate total magnitudes 
(valid for compact sources), and we have also corrected for 
the small amount of foreground Galactic extinction toward 
these fields using the COBE/DIRBE & IRAS/ISSA dust 
maps of Schlegel, Finkbeiner & Davis (1998). The opti- 
cal reddening is E(B — V) = 0.009, equivalent to extinc- 
tions of Assoip = 0.012, Aios™ = 0.010, A 125w = 0.008 & 
Ai 60 » = 0.005. 

2.3 Construction of Catalogues 

To perform the candidate selection we used the SExtrac- 
tor photometry package (Bertin & Arnouts 1996), version 
2.5.0. For Y-drops (objects clearly detected in the WFC3 J- 
band but with minimal flux in the Y-band and ACS images), 
apertures were 'trained' in the //-band image, and running 
SExtractor in dual-image mode those apertures were used 
to measure the flux in the same locations in the Y-band and 
J-band images. For each waveband we used a weight image 
derived from the exposure map. The z'-drop selection was 
done from catalogs trained in the J-band rather than in the 
//-band 

Tables [2] and [3] present our photometry of z'- and Y- 
drops from SExtractor. The MULTIDRIZZLE geometric 
transformation and image re-gridding produces an output 
where the noise is highly correlated, hence measuring the 
standard deviation in blank areas of the final drizzled image 
will underestimate the noise (e.g., Casertano et al. 2000). As 
in our previous work we have corrected the magnitude errors 
returned by SExtractor using our "true noise frames" , com- 
binations of the data obtained without using MULTIDRIZ- 
ZLE and hence without correlation between adjacent pixels, 
to determine the scaling factor (typically SExtractor un- 
derestimated the magnitude errors by a factor of ~ 2 for 
pixfrac=0.8 used for most of our data). We also measure 

1 Sec http://candels.ucolick.org/data_access/GODDS-S.html 
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WFC3 exposure times in ksec (5cr depth, AB mag) Area 
Field ID F-band J- band _ff-band (arcmin 2 ) 



CANDELS DEEP 8.1 (27.8) 7.4 (27.3) 7.7 (27.2) 62.9 

CANDELS WIDE 2.7 (26.8) 2.1 (26.9) 2.1 (26.6) 32.8 

Table 1. The total exposure time (in ksec) is listed for each WFC3 filter used in this study for both CANDELS 'wide' and 'deep' fields. 
In parenthesis, the average depth for each filter over the area listed is shown. These are 5<r limits calculated in apertures of CK'6 diameter, 
corrected as described in the text for aperture loss and reddening. 



the correlated noise (the standard deviation of the back- 
ground counts) in the drizzled image mosaics which we use 
for our source detection and photometry, and use the rela- 
tions in equation A13 of Casertano et al. (2000) to intro- 
duce a correction factor which depends on the output pixel 
scale and the size of the "droplet" in the drizzling proce- 
dure ("pixfrac"). We generally found good agreement (at 
the 0.05 mag level) with our sensitivity measurements using 
the true-noise frames. The errors displayed in Tables [5] and 
Table [3] are the corrected output from SExtractor 



3 CANDIDATE SELECTION 

Identification of candidates is achieved using the Lyman 
break technique (e.g., Steidel et al. 1996), where a large 
colour decrement is observed between filters either side of 
Lyman-a in the rest-frame of the galaxy. At z > 6, the flux 
decrement comes principally from the large integrated opti- 
cal depth of the intervening absorbers (the Lyman-a forest). 

At z ~ 8 — 9 the location of the Lyman-a break is 
redshifted to ~ 1.1/im - the WFC3 Yiosw and Ji25w are 
suitably located such that a 7.6 < z < 9.8 star forming 
galaxy will experience a significant flux decrement between 
these two filters, while for z ~ 7 the break lies at ~ 1/xm, 
between filters WFC3 Yiosw and ACS zssoip, with a redshift 
range of 6.5 < z < 8.0 (see Figure[T}. The selection efficiency 
drops at the extremes of these ranges. 

3.1 Selection Criteria 

Our photometrically selected Lyman-break sample suffers 
from contamination due to photometric scatter and interlop- 
ers (in particular L and T type dwarf stars and red galaxies 
at intermediate redshift). To discriminate candidates from 
these interlopers, we use the photometric data from another 
filter at wavelengths longer than the break, Ji25w for z'- 
drops and .Hi60w for K-drops and impose limitations on the 
2850i P — Vi05w and Ji25w — #i60w colours (respectively) as 
well, drawing a selection window in the colour-colour dia- 
gram that excludes most of the contaminants (Figure [5J). 

In this work we present objects within the colour - 
colour windows we selected and with detections of at least 
5 a in the two bands at wavelengths longer than the Lyman- 
a break. Even though the selection windows rule out most 
of the intrinsically red interlopers, these can still be included 
in our selection because of photometric scatter. To minimise 
this contamination, all objects with a > 2a detection in any 
of the 6435W, f606w and i775i p (below the Lyman limit) are 
classified as contaminants, ruling out in this way lower red- 
shift red galaxies, which we expected to faintly detect in the 
optical bands (see Figure [TJ. 



3.6 3.8 4 4.2 



50 Myr Corn. SF. Z=Z , z=8.0 




3.6 3.8 4 4.2 

log 10 (A[A]) (observed frame) 



Figure 1. Top panel - Model (from the Starburst99, Lcithcrcr 
et al. 1999) spectral energy distribution (SED) of a redshifted 
z = 8 star forming galaxy. Middle panel - Potential contaminants: 
Observed SED of a low-mass dwarf star (class: T4.5, Knapp et 
al. 2004) together with the model (Starburst99) SED of a 3.5Gyr 
Single-aged Stellar Population (SSP) at z = 2.5. The bottom 
two panels show the transmission functions of the combination of 
filters available to each field. 

Various colour selection windows have been proposed 
in the literature to remove contaminants and select high- 
redshift Lyman break galaxies. In this paper we use the cri- 
teria we derived previously for the z'-drops at z ~ 7 (Wilkins 
et al. 2011a, Wll hereafter): 

(Z850lp — ^105™) > 1-0 

{zs50i P - Y 105w ) > 2.4 x (Y 105w - Jus™) + 0.9 

(Vl05ui — Jl25u>) < 1-0 

and for the F-drops at z ~ 8 (Lorenzoni et al. 2011, Lll 
hereafter): 

(^105™ — Jl25uj) > 0.9 

{Y]_05 W — Jl25uj) > 0.73 X (Jl25iu — HifMw) + 0.9 
(Jl25iu — #160u>) < 1-5 

We also derive a list of candidates obeying the colour-cuts 
proposed by Bouwens et al. (2011) for these redshifts (we 
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la limit based on the noise and measured flux within the 
aperture. 




3.1.1 



-drops 



fl25w ''F160w 



In the 'deep' area we find 17 objects meeting our selection 
criteria (see Figure O top panel). Of these, 16 candidates 
meet the Bllz selection window, while 10 meet Wll (9 
of which also match the Bllz window). One of these ob- 
jects is UDFz-4256656 from Bouwens et al. (2011) in the 
HUDF field. In the 'wide' area, 2 objects meet the Bllz 
window (GS.W-zDl and GS.W-zD2), and none fall within 
the Wll colour selection. Images of the z'-drops meeting 
the selection criteria (Bllz and/or Wll) are shown in Fig- 
ure [3] Having three filters longwards of the break, it is 
possible to determine the UV spectral slope for the z - 
drop candidates: the Y-band filter could be affected by ei- 
ther the Lyman break or Lyman-a emission, or both, so 
the J- and //-bands are necessary to have 'clean' informa- 
tion on the UV slope. As in Wilkins et al. 2011b, /3 is de- 
termined from the (J125™ — Hieow) colour by the relation 
/3 = 4.28 x (Ji25u, — -ffieo™) — 2.0, which assumes that the 
slope is represented exactly by a power law. The /3 values 
are listed in Table [2] as already observed by Bunker et al. 
(2010), Wilkins et al. (2011b), Bouwens et al. (2010), the 
UV slopes of high-redshift galaxy candidates are very blue 
(/3 ~ —2), with fainter objects being bluer than the brighter. 
Note that the errorbars for faint candidates, due to photo- 
metric scatter, are considerable. 



3.1.2 Y -drops 

In the CANDELS 'deep' area, Table [3] 2 objects meet the 
Lll colour selection (Figure Q, both of which are included 
in the 6 objects selected with the BUY criteria (Figure [21 
bottom panel). We did not find any Y-drop candidate in 
CANDELS 'wide' area. 



Figure 2. Colour - colour diagrams for both z'-drops (top) and 
Y-drops (bottom). The shaded areas are the selection windows 
used, defined in Section 3.1 (light shading for the Bllz and BUY 
selection windows, darker shading for Wll and Lll). The objects 
we found are shown as grey dots (objects not meeting any of the 
colour - colour windows we are considering), black dots (objects 
in Bllz or BUY) and black circled dots (objects meeting Wll or 
Lll). The coloured dots denote the position of potential L and T 
dwarfs stars contaminants. The solid red line shows the colours 
that lower redshift galaxies (modelled as an instantaneous burst of 
star formation at z = 20 and no dust) would have, and the dotted, 
dashed and dot-dashed lines show this low-redshift template with 
reddenings of E(B -V) = 0.1, 0.25 & 0.5 respectively. The blue 
line is the predicted path taken by high-redshift galaxies (constant 
star formation from z = 20, no dust). For the red and blue tracks, 
numbers in correspondence with open circles indicate the redshift. 



label these Bouwens et al. criteria Bllz for z ~ 7 and BUY 
for z ~ 8 hereafter). This will allow for an easier comparison 
of candidates, and to investigate the effect of different selec- 
tion windows on the derivation of a luminosity function. For 
detections of less than la in the z' - or Y-band, we quote a 



3.1.3 Comparison to Other Studies 

Both the 'deep' and 'wide' CANDELS observations of 
GOODS-S have been recently searched for Y-drop candi- 
dates by both Oesch et al. (2012, hereafter 012) and Yan et 
al. (2012, hereafter Y12), resulting in 11 and 8 high redshift 
galaxy candidates respectively. 

Of the 11 012 sources we match only 3 with our 6 can- 
didates. Another object in our sample (GS.D-YD3) is also 
flagged as a potential candidate by 012 (CAND-2253348542) 
though is dismissed by 012 on the grounds of its stellar- like 
profile. We also match an additional 2 of our candidates with 
the 8 Y12 sources (there are no matches in common between 
all three candidate lists), thus all our candidates exist in ei- 
ther 012 or Y12. 

Given the lack of agreement between the previous catalogs 
of Y-drops (012 and Y12) with our new selection, and also 
the poor agreement between 012 and Y12 (there are 2 ob- 
jects in common of which neither is in our candidate list) it 
is useful to examine each of the 012 and Y12 candidates in 
turn to identify why they were not selected by us. 
Of the 8 012 sources not selected as candidates by us, 2 ob- 
jects, CANDY-2499448181 (which is also 048 in Y12) and 
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CANDY-2209651371, are detected at > 2a in a single optical 
band (though at < 3a). One object (CANDY-2320345371) 
is excluded because its (Y — J) colour is slightly bluer than 
our selection window, while a further 4 sources fail to meet 
our S/N> 5 criteria though do appear to be real objects (all 
detected at > 4a in both J- and ff-band). A single source 
(CANDY-219147298) is not matched within 0'.'5 of an ob- 
ject in our catalogue. Of the 6 Y12 sources not matched to 
our candidates the two brightest (048 and 100) are excluded 
on the basis of weak (2 to 3<r) optical detections in a single 
band. The 4 remaining objects are excluded on the basis of 
S/N concerns (in that they fall below S/N= 5 in one or both 
bands); in three cases (094, 035, 043) we detect the source 
at > 4cr in both Ji25w and -Hieow while the final object (085) 
is only detected at 2 — 3a and has colours inconsistent with 
our selection window. 

There are then two principal reasons for the Y12 and 012 
objects being excluded from our candidate list; at the bright- 
end 2 objects in each study (with 1 in common) are excluded 
due to weak (2 to 3a) optical detections in single band; while 
at the faint end several sources are excluded on the basis of 
our S/N criteria. In all but one case (Y12: 085) these objects 
are detected at > 4a in both Ji25w and -Hi60w and have ob- 
served colours consistent with our selection window. It then 
seems possible that some of the additional Y12 and 012 
candidates are potential high-redshift star forming galaxies. 
However, these objects are nevertheless excluded from the 
subsequent analysis of the rest-frame UV luminosity func- 
tion, as we want a robust sample. The computation of the 
effective volume takes into account our more conservative se- 
lection criteria, which should lead to the accurate luminosity 
function being recovered. 



4 DISCUSSION 

4.1 The bright-end of the UV luminosity function 
at z « 7 - 9 from CANDELS 

From our selection of z'- and Y-drops we can recover the vol- 
ume density of galaxies at z w 7 and z ~ 8 as a function of 
the rest-frame UV luminosity. The Lyman-break technique 
does not have uniform sensitivity on the probed redshift 
range, so we quantify the probability of recovering a high- 
redshift galaxy in our survey as a function of redshift and 
absolute UV magnitudes, p(Muv, z), with simulations. To 
perform these simulations we add into the images a large 
number of fake galaxies, with properties similar to those 
of the observed high-redshift population (i.e. compact with 
half-light radii rni ~ f/.'l, large Lyman-a forest decrement 
of Da ~ 0.99 and blue rest-frame UV colours). We then run 
our selection procedure and infer the probability of recover- 
ing such galaxies as a function of redshift and magnitude. 
From this probability the effective survey volume V e ff can 
be calculated, with the same approach described in Steidel 
et al. (1999) and Stanway, Bunker & McMahon (2003). We 
assume the LF to have a Schechter (1976) profile with four 
fixed values for a, -1.5, -1.7, -1.9 and -2.1, as the faint end 
slope cannot be strongly constrained with current data. The 
other Schechter parameters, </>* and M* 600 , are determined 
by maximising the Poissonian likelihood of observing a num- 
ber of objects in a magnitude bin. 




Figure 3. 2'' 4 X 2'.'4 bvizYJH thumbnail images of potential z at 
7 objects meeting our selection criteria in CANDELS GOODS- 
South field, ordered by J-band magnitude (brightest at the top). 
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z 7 



ID RA Doc z AB Y AB J AB H ab (z- (Y— p 

(J2000) (J2000) Y) AB J) AB 



GS 


D-zDl 


03 


32 


55 


930 


-27; 


19 


38, 


59 


26.98 ± 0.21 


26.01 ± 0.044 


25.87 ± 0.071 


26 


.25 


± 





.15 


0.97 


0.14 


—3.63 ± 0.74 


GS 


D-ZD2 1 


03 


:32 


37 


181 


-27: 


48: 


56 


68 


28.88 ± 0.74 


26.52 ± 0.043 


26.34 ± 0.074 


26 


46 


± 





12 


2.36 


0.18 


—2.51 ± 0.64 


GS 


D-zD3 1,2 


03 


:32 


08 


130 


-27: 


46: 


40, 


.88 


> 28.49 


26.85 ± 0.12 


26.37 ± 0.10 


26 


39 


± 





15 


> 1.64 


0.48 


—2.09 ± 0.82 


GS 


D-ZD4 1 


03 


32 


36 


006 


-27: 


11: 


11 


7 1 


> 28.26 


26.57 ± 0.067 


26.39 ± 0.085 


26 


.52 


± 





.14 


> 1.69 


0.18 


—2.56 ± 0.74 


GS 


D-zD5 


03 


32 


25 


447 


-27: 


50: 


53. 


36 


27.76 ± 0.33 


26.76 ± 0.094 


26.46 ± 0.084 


26 


11 


± 





.12 


1.0 


0.3 


— 1.91 ± 0.67 


GS 


D-zD6 


03 


32 


09 


583 


-27: 


46: 


32. 


06 


28.01 ± 0.42 


26 99 + 1 3 


26 64 + 11 


26 


01 


± 





.10 


1.02 


0.35 


70 + 73* 


GS 


W-zDl 


03 


32 


57 


.390 


-27: 


53: 


21 


77 


27.51 ± 0.25 


26.56 ± 0.18 


26.67 ± 0.16 


26 


,84 


± 





,24 


0.95 


—0.11 


—2.7 ± 1.2 


GS 


,W-zD2 


03 


:32 


36 


729 


-27: 


54: 


42. 


12 


27 35 + 21 


26 59 + 1 7 


26 77 + 14 


26 


58 


± 





17 


0.76 


—0.18 


— 1 19 + 94 


GS 


D-zD7 


03 


32 


36 


240 


-27: 


46: 


31 


37 


28.51 ± 0.65 


27.26 ± 0.12 


26.82 ± 0.11 


26 


78 


± 





.15 


1.25 


0.44 


—1.83 ± 0.85 


GS 


D-ZD8 1 


03 


:32 


10 


,693 


-27: 


41: 


16 


72 


> 28.09 


27.02 ±0.11 


26.83 ± 0.12 


26 


,71 


± 





16 


> 1.07 


0.19 


-1.49 ±0.91 


GS 


D-ZD9 1 


03 


:32 


:28 


.859 


-27: 


49: 


12 


63 


> 28.35 


27.18 ±0.12 


26.89 ±0.14 


26 


,88 


± 





20 


> 1.17 


0.29 


-1.96 ± 1.10 


GS 


D-zDlO 1 


03 


:32 


:27 


.916 


-27: 


45: 


42. 


72 


> 29.28 


27.24 ±0.17 


26.9 ±0.14 


27 


80 


± 





47 


> 2.04 


0.34 


-5.8 ±2.1* 


GS 


D-zDll 1 


03 


:32 


19 


,938 


-27: 


47: 


10. 


57 


29.01 ± 1.05 


27.04 ±0.10 


26.95 ± 0.13 


27 


,59 


± 





35 


1.97 


0.09 


-4.7 ± 1.6 


GS 


D-ZD12 1 


03 


:32 


17 


,638 


-27: 


48: 


29. 


21 


28.98 ± 0.93 


27.17 ±0.11 


27.07 ±0.15 


27 


66 


± 





39 


1.81 


0.1 


-4.5 ± 1.8 


GS 


D-zD13 


03 


:32 


12 


,512 


-27: 


47: 


56. 


86 


28.12 ±0.37 


27.22 ±0.12 


27.14 ±0.16 


27 


87 


± 





46 


0.9 


0.08 


-5.1 ± 2.1 


GS 


D-zDU 1 


03 


.32 


:37 


.230 


-27: 


45: 


38. 


41 


28.05 ± 0.44 


27.02 ±0.10 


27.15 ± 0.15 


27 


36 


± 





26 


1.03 


-0.13 


-2.9 ± 1.3 


GS 


D-ZD15 1 


03 


:32 


:30 


.793 


-27: 


50: 


27 


19 


> 29.03 


27.46 ±0.17 


27.17 ±0.15 


27 


56 


± 





32 


> 1.57 


0.29 


-3.7± 1.6 


GS 


D-zD16 


03 


:32 


16 


.057 


-27: 


47: 


57. 


72 


28.09 ± 0.44 


27.3 ±0.13 


27.21 ± 0.16 


27 


.70 


± 





37 


0.79 


0.09 


-4.1 ± 1.8 


GS 


D-zD17 


03 


:32 


:35 


.067 


-27: 


46: 


34 


96 


28.35 ±0.52 


27.51 ±0.15 


27.22 ± 0.15 


27 


89 


± 





11 


0.84 


0.29 


-4.8 ±2.0* 



in Wll selection.; 2 not selected using Bll criteria.; * outside the colour-colour selection window employed by Wilkins et al. (2011b) 

for a clean selection of 2-drops for analysis of spectral slope, p. 



Table 2. z'-band drop out candidate at z f» 7 meeting either of the selection criteria described. Objects are ordered by apparent Jab 
magnitude. Where quoted, limits are la 



z 8 



ID 


RA 


Dec 


Y AB 


Jab 


H AB 


(Y - J) AB 


{J-H) AB 


Bll 


Lll 


GS.D-YD1 


03:32:48.921 


-27:47:07.36 


27.0 ±0.11 


26.18 ±0.063 


26.17 ±0.077 


0.82 


0.01 


/ 




GS.D-YD2 


03:32:14.135 


-27:48:28.96 


28.18 ±0.3 


26.94 ±0.12 


26.8 ±0.13 


1.24 


0.14 


/ 


/ 


GS.D-YD3 


03:32:25.330 


-27:48:54.07 


27.18 ±0.11 


26.59 ±0.086 


26.9 ± 0.13 


0.59 


-0.31 


/ 




GS.D-YD4 


03:32:44.018 


-27:47:27.23 


27.8 ±0.19 


27.01 ± 0.13 


26.97 ±0.16 


0.79 


0.04 


/ 




GS.D-YD5 


03:32:40.257 


-27:44:09.84 


27.61 ±0.18 


27.09 ±0.14 


27.02 ±0.16 


0.52 


0.07 


/ 




GS.D-YD6 


03:32:20.979 


-27:48:53.46 


29.04 ±0.64 


27.0 ±0.13 


27.05 ±0.16 


2.04 


-0.05 


/ 


/ 



Table 3. Y-band drop out candidate at z m 8 meeting either of the selection criteria described. Objects are ordered by apparent H AB 
magnitude. 



In Figures [5] and [6] we plot our datapoints at z ~ 7 
(Wll selection window) and z « 8 (Lll selection window), 
respectively, against several luminosity functions from our 
previous work (Wilkins et al. 2011a, Lorenzoni et al. 2011) 
and other publications (Bouwens et al. 2011, Oesch et al. 
2012, Yan et al. 2010). In the same Figures we also plot our 
datapoints obtained for the Bllz and BUY selection win- 
dows. As can be clearly seen, the number densities inferred 
from the different selection windows are in good agreement, 
within the error bars. We will therefore consider the lists 
of candidates obtained using the Bllz and BUY selection 
windows. 

In Tableland Table [6] we show the best fitting results 
for the LF at redshifts z ~ 7 and z ~ 8, respectively, for each 
of the selections windows used. The candidates found in our 
previous works in the HUDF and ERS fields (Wilkins et 
al. 2011a, Lorenzoni et al. 2011) are also included in all the 
LF calculations. In fitting the Schechter luminosity function, 
phi* and M* are highly correlated, so we show the error 



ellipses (1 a and 2 a significance contours) for the z ~ 7 and 
2 ~ 8 — 9 luminosity functions in Figure [7] 

We note very good agreement at z ~ 7 between the 
best fitting LFs obtained using the two different selection 
windows (Wll and Bllz). These results are also in line with 
our previous estimates (Wilkins et al. 2011a). 

At z w 8, the Lll selection window adds only 2 can- 
didates to our previous sample of y-drops. The BUY se- 
lection yields 6 candidates, and combining these with our 
previous sample (accounting for the different effective vol- 
umes probed by the colour selections) produces luminosity 
functions (Table [B]) consistent with those of several previ- 
ous studies Lorenzoni et al. (2011), Bouwens et al. (2011), 
Oesch et al. (2012), which indicate fainter characteristic lu- 
minosity, L* , than at lower redshifts. However, these results 
at z ~ 8 are strongly inconsistent with the LF proposed by 
Yan et al. (2010) on the basis of their analysis of the HUDF, 
in which they claimed far more faint y~drop galaxies than 
in the analyses of other groups (Bunker et al. 2010; McLure 
et al. 2010; Bouwens at el. 2010). As can be seen in Figure[6l 



© 0000 RAS, MNRAS 000, 000-000 




w ; 



.D-YD3 



The Bright-end of the UV LF at z « 7 - 9 /rom CANDELS 7 













rr HE - . .-. I I 





Figure 4. 2" 4 X 2'.'4 bvizYJH thumbnail images of potential z Rj 
8 objects meeting our selection criteria in CANDELS GOODS- 
South field, ordered by H-band magnitude (brightest at the top). 
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Yan '12 Class 



GS.D-YD1 
GS.D-YD2 
GS.D-YD3 
GS.D-YD4 
GS.D-YD5 
GS.D-YD6 



CANDY- 
CANDY 
CANDY 

CANDY 
CANDY 
CANDY 
CANDY 
CANDY 
CANDY 
CANDY- 
CANDY- 
CANDY- 



2141348289 
2253348542 
2440247273 

2209848535 
2499448181 
2320345371 
2209651371 
2350049216 
2192147298 
2181852456 
2379552208 
2408551569 



064 

107 
048 

035 



100 

094 
043 
085 



O 

w 
o 

F 

? 

F 
F 
F 
O 
F 
F 
F, W 



F - Object too faint in J— and/or H— band for our selec- 
tion criteria. 

O - Detection of more than 2 a in at least one of the optical 
bands. 

W - Object outside our colour-colour selection windows. 
? - Object not picked up by SExtractor. 



Table 4. We list here candidates identified by Oesch et al. (2012) 
and Yan et al. (2012), second and third column respectively, and 
match them with ours when possible (first column) or give the 
reason why we do not find them (fourth column). 



our measured number densities of Y-drops at brighter mag- 
nitudes (Muv = —21 & —20) are inconsistent by an order 
of magnitude or more than the expectation from the Yan et 
al. (2010) LF. 

We now compare the star formation rate (SFR) densi- 
ties obtained by integrating the z ~ 7 and z « 8 luminosity 
functions down to various limiting magnitudes (Figures [8] 



Figure 5. The luminosity distribution (top) and luminosity func- 
tion (bottom) of z'-drop selected sources at z a 7. Our datapoints 
are plotted against Wilkins et al. (2011a, solid line) and Bouwcns 
et al. (2011, dashed line) luminosity functions. The uncertainty 
bars represent the 68.2% poissonian confidence interval of the 
number density <j>. The upper limits denote the maximum value 
of the 68.2% confidence interval with n = observations. This 
corresponds roughly to n = 1.84, i.e. the for an observed n = 
there is a 68.2% chance the true value is < 1.84. 



&[9]to the SFR densities required for reionization from the 
Madau, Haardt & Rees (1999) relation: 



PSFR 



0.012 M Q yr" 1 Mpc" 

/esc 



1 + 2 

1 + 8.6 



0.0462 



We have updated equation 27 of Madau, Haardt & Rees 
(1999) for a more recent concordance cosmology estimate 
of the baryon density from Larson et al. (2010), €lb /iioo = 
0.022622. In the above equation, C is the clumping factor 
of neutral hydrogen, C = (pm) {pm)~ 2 , whose used value 
in this work is 5 (Pawlik et al. 2009). / OS c is the escape 
fraction of ionizing photons, which is highly uncertain - we 
consider escape fractions as high as 100 per cent (rather 
implausible) and down to 10 per cent (which may be the 
average at z ~ 3 population, Nestor et al. 2011). At z ~ 8.6 
(the average redshift of the y-drops) , reionization cannot be 
achieved with the observed luminosity functions unless the 
slope is a = —1.9 or steeper, even if the escape fraction is 100 
per cent. However, a steeper faint end slope, an even lower 
IGM clumping factor, and a low-metallicity population (or 
a top-heavy IMF) might still provide sufficient photons for 
star-forming galaxies to reionize the Universe (see Lorenzoni 
et al. 2011). 
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z 7 

Wll Bllz 



a 


M* 600 [AB mag] 


<j>* [Mpc- 3 ] 


M 1600 I AB ma g] 


[Mpc~ 3 ] 


— 1.5 


-19.75 


0.00152 


-19.75 


0.00159 


— 1.7 


-19.95 


0.00110 


-19.93 


0.00119 


-1.9 


-20.19 


0.00072 


-20.14 


0.00081 


-2.1 


-20.51 


0.00039 


-20.40 


0.00049 



Table 5. The best fit values for M^qq and <f>* at z ss 7 for a Schechter function assuming fixed a £ {—1.5, —1.7, —1.9, —2.1} for both 
the Wll (columns 2 and 3) and Bllz (columns 4 and 5) selection windows. 



z 8 



Lll BUY 



Q 


M l*600 I AB ma §] 


</>* [Mpc~ 3 ] 


M l*600 I AB ma §] 


[Mpc~ 3 ] 


— 1.5 


-19.10 


0.00143 


-19.42 


0.00088 


— 1.7 


-19.23 


0.00119 


-19.53 


0.00075 


-1.9 


-19.37 


0.00095 


-19.66 


0.00060 


-2.1 


-19.54 


0.00069 


-19.80 


0.00046 



Table 6. The best fit values for Afj*gQ and (f>* at z fa 8 for a Schechter function assuming fixed a £ {—1.5, —1.7, —1.9, —2.1} for both 
the Lll (columns 2 and 3) and BUY (columns 4 and 5) selection windows. 
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Figure 7. The likelihood contours for the luminosity function of 
,z'-drops (dashed lines) and Y-drops (solid lines), showing the cor- 
relation between the fitted M* and <j>* parameters for a Schechter 
function fit, using our sample of galaxies from the Bll colour se- 
lection. A faint-end slope of a = —1.9 is adopted here. The 68% 
(inner) and 95% (outer) likelihood contours are shown.The cross 
represents the best-fit parameter values. 



Figure 6. The luminosity distribution (top) and luminosity func- 
tion (bottom) of Y-drop selected sources at z xs 8. Our data- 
points arc plotted against several luminosity functions: Lorenzoni 
et al. (2011, solid dark line), Oesch et al. (2012, solid light line), 
Bouwens et al. (2011, dashed line) and Yan et al. (2010, dotted 
line). The uncertainty bars represent the 68.2% poissonian confi- 
dence interval of the number density <j>. The upper limits denote 
the maximum value of the 68.2% confidence interval with n = 
observations. This corresponds roughly to n = 1.84, i.e. the for an 
observed n = there is a 68.2% chance the true value is < 1.84. 



5 CONCLUSIONS 

In this paper we present a list of candidate high-redshift star- 
forming galaxies identified with the Lyman-break technique 
using HST/WFC3 near-infrared data within the CANDELS 
programme. We have presented the first analysis of z'-drop 
candidate galaxies at z ~ 7 images with HST/WFC3 in the 
new CANDELS imaging of the GOODS-S field, building on 
previous work by our team (Bunker et al. 2010; Wilkins 
et al. 2010, 2011a) in the smaller HUDF and ERS fields 
within GOODS-S. We also use the colour selections derived 
by Lorenzoni et al. (2011) and Bouwens et al. (2011) to 
identify candidate z w 8 Y-drops galaxies in this field, and 
compare our catalogues with those independently derived 
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- 100 
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Figure 8. The solid lines are the total star formation rate den- 
sity (left axis) or ionising flux density (right axis) inferred from 
the luminosity function fits for our x'-drop sample (for faint end 
slopes a = [—1.5, —1.7, —1.9, —2.1]), integrating down to the lim- 
iting absolute magnitude in the rest-frame UV shown on the lower 
x-axis (in AB magnitudes); the upper x-axis shows the equiva- 
lent unobscured star formation rate. The dashed lines show the 
requirement to keep the Universe ionised at z — 7, using the re- 
lation from Madau, Haardt & Rees (1999) and assuming a low 
clumping factor of C = 5. We show the requirements for escape 
fractions of / esc = 0.1, 0.5 & 1. Where the solid lines cross the 
dashed lines, reionzation can be achieved. The shaded region is 
where the current deepest observations probe (the HUDF). 
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Figure 9. The solid lines are the total star formation rate den- 
sity (left axis) or ionising flux density (right axis) inferred from 
the luminosity function fits for our Y-drop sample (for faint end 
slopes a = [—1.5, —1.7, —1.9, —2.1]), integrating down to the lim- 
iting absolute magnitude in the rest-frame UV shown on the lower 
x-axis (in AB magnitudes); the upper x-axis shows the equiva- 
lent unobscured star formation rate. The dashed lines show the 
requirement to keep the Universe ionised at z = 8.6, using the 
relation from Madau, Haardt & Rees (1999) and assuming a low 
clumping factor of C = 5. We show the requirements for escape 
fractions of / csc = 0.1, 0.5 & 1. Where the solid lines cross the 
dashed lines, reionzation can be achieved. The shaded region is 
where the current deepest observations probe (the HUDF). 



from the same CANDELS field by Oesch et al. (2012) and 
Yan et al. (2012). We treble the number of bright (H mag < 
27) Y-drops from Lorenzoni et al. (2011) and double the 
number of bright (J mag < 27.2) z'-drops from Wilkins et al. 
, (2011a). 

The bright high redshift galaxy candidates we found 
serve to better constrain the bright end of the luminosity 
function at those redshift, and may also be more amenable to 
spectroscopic confirmation than the fainter ones presented 
in various previous work on the smaller fields (HUDF and 
ERS). Indeed, with AB magnitudes of ~ 26 (longward of 
the break) , we could hope to detect Lyman-a emission lines 
with rest frame equivalent widths of a few tens of Angstroms 
(typical of Lyman break galaxies at z ~ 3 — 6, e.g. Stanway 
et al. 2004) in ~ 5 hours' spectroscopy with an instrument 
such as XSHOOTER on VLT (see Caruana et al. 2012). If 
spectroscopy reveals that Lyman-a does not emerge at these 
redshifts, then our bright Lyman-break galaxy sample can 
potentially place strong constraints on the absorption of the 
Gunn-Peterson (1965) damping wing (and hence the neutral 
fraction of hydrogen at z ~ 8). 

We also look at the agreement with previous luminos- 
ity functions derived from WFC3 drop-out counts, and find 
good agreement with those of Wilkins et al. (2011a) and 
Bouwens et al. (2011) at z ~ 7, and Lorenzoni et al. (2011) 
and Oesch et al. (2012) at z ~ 8. However, our results 
strongly rule out the z ~ 8 luminosity function proposed 
by Yan et al. (2010). 
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